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ABSTRACT 

^ ■ We have assembled a large sample of virialized systems, comprising 66 galaxy clusters, 

' groups and elliptical galaxies with high quality X-ray data. To each system we have 

, fitted analytical profiles describing the gas density and temperature variation with 

' radius, corrected for the effects of central gas cooling. We present an analysis of the 

scaling properties of these systems and focus in this paper on the gas distribution 
CD , and M — Tx relation. In addition to clusters and groups, our sample includes two 

early-type galaxies, carefully selected to avoid contamination from group or cluster X- 
ray emission. We compare the properties of these objects with those of more massive 
systems and find evidence for a systematic difference between galaxy-sized haloes and 
Qh' groups of a similar temperature. 

JL , We derive a mean logarithmic slope of the M — Tx relation within i?200 of 

5_, ■ 1.84±0.06, although there is some evidence of a gradual steepening in the M — Tx re- 

I lation, with decreasing mass. We recover a similar slope using two additional methods 

, of calculating the mean temperature. Repeating the analysis with the assumption of 

isothermality, we find the slope changes only slightly, to 1.89±0.04, but the normaliza- 
tion is increased by 30 per cent. Correspondingly, the mean gas fraction within i?200 

_ 3 _ 3 

changes from (0.13 ± 0.01)h^if to (0.11 ± 0.01)h^if , for the isothermal case, with the 
smaller fractional change reflecting different behaviour between hot and cool systems. 
There is a strong correlation between the gas fraction within 0.3i?2oo and temperature. 
This reflects the strong (S.Scr) trend between the gas density slope parameter, (3, and 
temperature, which has been found in previous work. 

These findings are interpreted as evidence for self-similarity breaking from galaxy 
feedback processes, AGN heating or possibly gas cooling. We discuss the implications 
of our results in the context of a hierarchical structure formation scenario. 

Key words: galaxies: clusters: general - galaxies: haloes - intergalactic medium - 
X-rays: galaxies - X-rays: galaxies: clusters 
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1 INTRODUCTION 

The formation of structure in the Universe is sensitive to 
physical processes which can influence the distribution of 
baryonic material, as well as cosmological factors which ulti- 
mately govern the behaviour of the underlying gravitational 
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potential. By studying the properties of groups and clus- 
ters of galaxies, it is possible to probe the physical processes 
which shape the evolution and growth of virialized systems. 

X-ray observations of the gaseous intergalactic medium 
(IGM) within a virialized system provide an ideal probe of 
the structure of the halo, since the gas smoothly traces the 
underlying gravitational potential. However, this material 
is also sensitive to the influence of physical processes aris- 
ing from the interactions between and within haloes, which 
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are commonplace in a hierarchically evolving universe (e.g. 
iBlumenthal et al1ll984ll . Even in relatively undisturbed sys- 
tems, feedback from the galaxy members can bias the gas 
distribution with respect to the dark matter in a way which 
vari es systematically with halo m ass. N-body simulations 
fe.g. lNavarro. Frenk fc Whitell995ll indicate that, in the ab- 
sence of such feedback mechanisms, the properties of the 
gas and dark matter in virialized haloes should scale self- 
similarly, except for a modest variation in dark matt er con- 
centration with mass jNavarro. Frenk fc Whitlll997l) . Con- 
sequently, observations of a departure from this simple ex- 
pectation provide a key tool for investigating the effects 
of non-gravitational heating mechanisms, arising from feed- 
back processes. 

There is now clear evidence that the properties of clus- 
ters and groups of galaxies do not scale self-similar ly: for 
example, the L — Tx relation in clusters sh ows a logarithmic 
slope which is steeper than expected (e.g. [ Edge fc Steward 
ll99ll:lArnaud fc Evrardll999l:lFairlev et all2000t) " A further 
steepening of this slope is observed in the group regime (e.g. 
iHelsdon fc: Po nman"2000V consistent with a flattening in 
the gas densit y profiles, which is evident in syste ms cooler 
than 3-4 keV JPonman. Cannon fc Navarroll999^ . Such be- 
haviour is attributed to the effects of non-gravitational heat- 
ing, which exert a disproportionately large influence on 
the smallest haloes. An obvious candidate for the source 
of this heating is galaxy winds, since these are known to 
be responsible for the enrichment of the IGM with heavy 
elements (e.g. iFinoeuenov. Arnaud fc DavidI lioOlf) . How- 
ever, active galactic nuclei (AGN) may also play a signif- 
icant role, particularly as there is some debate over the 
amount of energy available from supernova-driven outflows 
iWu. Fabian fc Nulsenl2'o0ol) . Recently, t heoretical work has 
also examined the role of gas cooling fc.f. lKnight fc PonmanI 
I1997I) . which is also able to reproduce the observed scaling 
properties of groups and clusters, by eliminating the low- 
est entropy gas throug h star formation, thus allowing hotter 
material to replace it jMuanwong et al ] |200ll:IVoit fc BrvanI 
l200lD . 

Previous observational studies of the distribution of 
matter within clusters hav e typically been limited by e i- 
ther a small sample size (e.g. lDavid. Jones fc^brman| l99^). 
or ha ve assumed an isothermal IGM (e.g. Iwhit^^^ibianl 
[W9^; it appears tha t significant tem perature gradients are 
present in many (e.g.|MaAevitcI^^^.,^^1998 ^'). although per- 
haps not aU fe.g. I WhitdbOOd : iDe Grandi fc Molendill2002t 
llrwin fc Bregmanl200(ir ' clusters of galaxies. Another issue is 
the restriction imposed by the arbitrary limits of the X-ray 
data; halo properties must be evaluated at constant frac- 
tions of the virial radius (-Rv), rather than at fixed metric 
radii imposed by the data limits, in order to make a fair 
comparison between varying mass scales. In this work, we 
derive analytical expressions for the gas density and temper- 
ature variation, which allow us to extrapolate these quanti- 
ties beyond the limits of the data. However, we are careful to 
consider the potential systematic bias associated with this 
process. Our study combines the benefits of a large sample 
with the advantages of a 3-dimensional, deprojection anal- 
ysis, in order to investigate the scaling properties of virial- 
ized haloes, spanning a wide range of masses. In this work 
we have brought together data from three large samples, 
comprising the majority of the suitable, radially-resolved 3D 



temperature analyses of clusters. We include a large num- 
ber of cool groups in our analysis, as the departure from 
self-similarity is most pronounced in haloes of this size: the 
non-gravitationally heated IGM is only weakly captured in 
the shallower potentials wells of these objects. 

To further extend the mass range of our analysis, we 
include two galaxy-sized haloes in our sample, in the form 
of an elliptical and an SO galaxy. Galaxy-sized haloes are 
of great interest as they represent the smallest mass scale 
for virialized systems and constitute the building blocks 
in a hierarchically evolving universe. Great emphasis was 
placed on identifying galaxies free of contamination from X- 
ray emission associated with a group or cluster potential, in 
which they may reside, since this is known to complicate 
analysis of their haloes (e.g. iMulchaev fc Zabludofll [19981 : 
iHelsdon fc PonmanI I2OO0I1 . The most well-studied galaxies 
are generally the first-ranked members in groups or clusters, 
and it is known that such objects are atypical, as a conse- 
quence of t he dense gaseous envir onment surrounding them: 
the work of IHelsdon et alJ (1200 ll) has shown that brightest- 
group galaxies exhibit properties which correlate with those 
of the group as a whole, possibly because many of them lie 
at the focus of a group cooling flow. The study of lSato et alJ 
i2000ll incorporated three ellipticals, but any X-ray emis- 
sion associated with these objects is clearly contaminated 
by emission from the group or cluster halo in which they are 
embedded. 

Throughout this paper we adopt the following cosmo- 
logical parameters; Ho = 70 km s~^ Mpc~^ and qo = 0. Un- 
less otherwise stated, all quoted errors are la on one param- 
eter. 



2 THE SAMPLE 

In order to investigate the scaling properties of virialized 
systems, we have chosen a sample which includes rich clus- 
ters, poorer clusters, groups and also two early-type galaxies, 
comprising 66 objects in total. Sample selection was based 
on two criteria: firstly, that a 3-dimensional gas temperature 
profile was available. In conjunction with the corresponding 
gas density profile, this allows the gravitating mass distribu- 
tion to be inferred. Secondly, we reject those systems with 
obvious evidence of substructure, where the assumption of 
hydrostatic equilibrium is not reasonable; it is known that 
the properties of such systems differ systematically from 
those of relaxed clusters (e.g. Ritc hie fc Tho mas 2002). This 
also favours the assumption of a spherically symmetric gas 
distribution, which is implicit in our deprojection analysis. 
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1.58 


M 


P,G,S 


re 


Abell 780 


0.0565 


4.63ts;ii 


2032l}g 


10.09l°:^2^ 


1-681H1 


0.67l°:°l 


- 


0, 


qr)+o.03 

■y^-0.03 


0.45 


L 


P+G 


Abell 2256 


0.0581 


8.62t»;^? 


18141}^^ 


q 1 + 0. 03 

'-'■-'^"-o.os 


5.02^rn 


0.78^ni 


- 


1, 


97+0.07 
■^'-0.05 


1.53 


M 


P,G,S 





Abell 1795 


0.0622 


8.54ti;«« 


ZOOOlf/o 


4.301^5 


4.011°:^? 


0.83^:22 


- 


1, 


1 7+0.10 
' -0.14 


1.44 


M 


P,G,S 




Abell 3112 


0.0703 


7.76lJ;«^ 


131lii5 


14.82lg;«^ 


l.UO_Q gg 


fis+° °^ 

U-O'J_0.02 


- 


1, 


09+0.14 
•"^^-0.09 


1.29 


M 


P,G,S 


re 


Abell 644 


0.0711 


11.68li:g 




7.76l°;g 


1 0+0. 18 


n 70+0.02 
' '-'-0.02 


- 


1, 


qr+O.ll 
•"^^-0.10 


1.27 


M 


P,G,S 




Abell 399 


0.0722 


7 07+0.69 
' ' -0.73 


1734li« 




1 sq+0-36 


A 1:0+0.05 

"•JJ-O.OS 


- 


1, 


lfi+009 
•^"-0.06 


1.26 


M 


P,G,S 




Abell 401 


0.0739 


9.55l°:f 


185lli^^3 




r\ q.7+0.09 

^•■J' -0.09 


n fiq+o oi 


- 


1, 


99+0.05 
•^^-0.04 


1.23 


M 


P,G,S 




Abell 2670 


0.0759 


5.64t°;^9 


16471 


6.20l°;l^ 


n 07+0 04 


A i:i;+0.02 
U.OO_Q 02 


- 


1, 


04+004 


- 


F 


P,S 




Abell 2029 


0.0766 




226611JJ 


fi Q4+0.10 
D-J^-O.lO 


r\ qy+0.09 

^"^'-0.09 


0.68^n3 


n on+0 06 




- 


1.69 


F 


P,S 




Abell 1650 


0.0845 




18161^?^ 


c c-l +0.55 
O-Oi_0.55 


9 oc^+0.78 
^•^^-0.78 


0.78^°:!^ 


- 


1, 


-|q+o.i5 

■^•^-0.17 


1.09 


M 


I,G,S 




Abell 1651 


0.0846 




1777lil° 


c oc+0.41 
"■^^-0.40 


9 ri9+0.23 
^•'-'^-0.23 


0.701^2 


- 


1, 


1 A+0.04 
•^^-0.05 


1.09 


M 


P,G,S 




Abell 2597 


0.0852 


r. Q^ + 0.47 


184lli«l 


6.471°;^^ 




0.68l°:g3 


— 


1 


Ar + 0.04 
•"-"^-0.04 


1.55 


F 


P,S 




Abell 478 


0.0882 




17231^3^ 


6.98l°:^l 


9 0/1+0.23 
^•■^^-0.23 


n 7c;+0-0i 


— 


1, 


04+0.17 
"J^-0.18 


1.06 


M 


P,G,S 




Abell 2142 


0.0894 




22161^2^^ 


r 91 +0.13 


14+0.22 
•J-J^^-0.22 


n 74+0.01 


— 


1, 


1 Q+0.10 
-L°-0.13 


1.05 


M 


P,G,S 




Abell 2218 


0.1710 


°-^°-1.33 


1 r\r\ ^-1-180 

1904+i°^ 


^ r7-(-0 15 


0.90* 


0.59l°:gi 




1 


■^J^-0.02 




L 


P+G 




Abell 665 


0.1818 


8.60ljjl 


r)t-jyo+268 


6 S^+O-i'^ 
"■■J^-o.ie 


1 91+0.04 
^•^-"^-0.04 


A fiK+O.Ol 




1, 


A9+O.O8 
■"^-0.06 




L 


P+G 




Abell 1689 


0.1840 


lz.o±_Q 93 


2955+i?| 


33.6li;;«* 


"•"'-'-0.00 


A 70+0.06 
"• ' 'J-0.00 


0.00* 






2.40 


L 


P+G 




Abell 2163 


0.2080 




2104l^«^ 


Y Y7+0.53 
'■' '-0.52 


1 63+°°** 


A 70+0.02 
"• ' 'J-D.02 




1, 


00+0.11 
"J"-0.22 


0.54 


M 


P,G,S 





Table 1: Som e key properties of the 66 objects in the sample, listed in order of increasing redshift. Redshifts are taken from 
lEbehng et all 11996. .1998,' l: .Ponman et alJ (11999) and NED. Columns 3-9 are data as determined in this work. All errors are 68% 
confidence. 

I indicates the two galaxies. 

* denotes no errors available, as parameter poorly constrained. 

also known as NGC 383. 
^ also known as Abell 4038. 

" The cooling-flow corrected, emission- weighted temperature of the system within 0.3ii2oo, as determined in this work. 
^ Temperature gradient; positive values mean T decreases with radius. 

" Cooling flow excision radius (M sample) or radius within which a cooling flow component was fltted (F,L,S samples) 

'' F = Finoguenov et al. , L = Lloyd-Davies et al. , M = Markevitch et al. , S = Sanderson et al. (this work) 

" P = ROSAT PSPC, H = ROSAT HRI, G = ASCA GIS, S = ASCA SIS, I = Emstem IPC; + denotes simuhaneous flt 
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By combining three samples from the work of Marke- 
vitch, Finoguenov and Lloyd-Davies (described in detail 
in sections 13.31 13.41 & 13.51 respectively) together with new 
analysis of an additional six targets (also described in sec- 
tion I^^J, we have assembled a large number of virialized 
objects with high-quality X-ray data. From these data, we 
have derived deprojected gas density and temperature pro- 
files for each object, thus freeing our analysis from the sim- 
plistic assumption of isothermality which is often used in 
studies of this nature. The large size of our sample ensures 
a good coverage of the wide range of emission- weighted gas 
temperatures, spanning 0.5 to 17 keV. Thus, we incorpo- 
rate the full range of sizes for virialized systems, down to 
the scale of individual galaxy haloes. The redshift range is 
z = 0.0036-0.208 (0.035 median), with only four targets ex- 
ceeding a redshift of 0.1. Some basic properties of the sample 
are summarised in Tabled 

As a number of systems are common to two or more of 
the sub-samples, we are able to directly compare data from 
different analyses, allowing us to investigate any systematic 
differences between the techniques employed. We present the 
results of these consistency checks in section 2] The diverse 
nature of our sample, with respect to the different methods 
used to determine the gas temperature and density profiles, 
insulates our study to an extent from the bieis caused by re- 
lying on a single approach. However, we are still able to treat 
the data in a homogeneous fashion, given the self-consistent 
manner in which the cluster models are parametrized (see 
section 13.11 . 



3 X-RAY DATA ANALYSIS 

The X-ray data used in this study were taken with the 
ROSAT PSPC and ASCA GIS & SIS instruments. Ahhough 
now superseded by the Chandra and XMM-Newton obser- 
vatories, these telescopes have extensive, publicly available 
data archives and are generally well-calibrated. In addition, 
the PSPC and GIS detectors have a wide field of view, which 
is essential for tracing X-ray emission out to large radii, par- 
ticularly for nearby systems, whose virial radii can exceed 
one degree on the sky. The use of three separate detectors, on 
two different telescopes, enhances the robustness of our anal- 
ysis, by reducing potential bias associated with instrument- 
related systematic effects. 

Since this work brings together data from separate sam- 
ples, there is considerable variation in the form in which 
those data were originally obtained. This necessitated a sup- 
plementary processing stage to convert the data into a uni- 
fied format, in order to treat them in a homogeneous fash- 
ion. In the case of the Finoguenov sample, analytical pro- 
files were fitted to deprojected gas density and temperature 
points (see section l^^ for details); for the Markevitch sample 
it was necessary to calculate the gas density normalization 
for such an analytical function, from the fitted data (sec- 
tion However, our chosen model parametrization - de- 
scribed below - was fitted directly to the raw X-ray data for 
the remaining systems, including the Lloyd-Davies sample 
(further details of the data analysis are given in section IX^ . 



3.1 Cluster models 

In order to evaluate the gas temperature and density in 
a virialized system, as well as derived quantities such 
as gravitating mass, at arbitrary radii, we require a 3- 
dimensional analytical description of these data. A core in- 
dex parametrization of the gas density, p(r), is used, such 
that 



p(r) = p(0) 



l+(^ 

Tc 



(1) 



where Tc and (5 are the density core radius and index pa- 
rameter, respectively. The motivation for the use of this 
parametrization is essentially empirical, although simula- 
tions of cluster mergers are capable of reproducing a core in 
the gas density, despite the cuspy nature of the underlying: 
dark matter distribution (e.g. lPearce. Thomas fc CouchmanI 
11994*). However, in the absence of merging, N-body simula- 
tions offer no clear explanation for the presence of a signifi- 
cant core in the IGM profile, even when t he effects of galaxy 
feedback mechanisms are incorporated (|Metzler fc EvrardI 
E997). 

The density profile is combined with an equivalent ex- 
pression for the temperature spatial variation, described by 
one of two models; a linear ramp, which is independent of 
the density profile, of the form 



T(r) = T(0) 



(2) 



where a is the temperature gradient. Alternatively, the tem- 
perature can be linked to the gas density, via a polytropic 
equation of state, which leads to 



T{t) = T(0) 



1 + 



1/3(7-1) 



(3) 



where 7 is the polytropic index and and /3 are as defined 
previously. 

Together, p(r) and T(r) can be used to determine the 
cluster gravitating mass profile as, in hydrostatic equilib- 
rium, the following condition is satisfied 



Mgrav (r) 



kT{r)r 



d In p d In T 



d In r d In r 



(4) 



^arazirJIlflSST ). where fi is the mean molecular weight of the 
gas and nip is the proton mass. This assumes a spherically 
symmetric mass distribution, which has been shown to be 
a reason able approximation, even for moderate ly elliptical 
systems jFabricant. Rvbicki fc Gorensteinlll984^ . 

Since the X-ray emissivity depends on the product of 
the electron and ion number densities, we parametrize the 
gas density in terms of a central electron number density 
(i.e. at r = 0), assuming a ratio of electrons to ions of 1.17. 
We base our inferred electron densities on the X-ray fiux 
normalized to the ROSAT PSPC instrument, as there is a 
known efi'ective area offset between this detector and the 
ASCA SIS and GIS instruments. In those systems where 
the original density normalization was defined differently, a 
conversion was necessary and this is described below. 

Once the gravitating mass profile is known (from equa- 
tion the corresponding density profile can be found triv- 
ially, given the spherical symmetry of the cluster models. 
This can then be converted to an overdensity profile, S{r), 
given by 
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5{r) = 



Ptot{r) 

pcrit 



(5) 



where ptotir) is the mean total density within a radius, r, 
and Pcrit is the critical density of the Universe, given by 
3Hi/8nG. 

It is the overdensity profile which determines the virial 
radius (-Rv) of the cluster; simulations indicate that a rea- 
sonable approxi mation to Ry is give n by the value of r when 
5{r) = 200 fe.g. lNavarro et al.|[l995ll - albeit for ptot{r) cal- 
culated at the redshift of formation, Z{, rather than the red- 
shift of observation, Zobs- and we adopt this definition in 
this work. Strictly speaking, the approximation = i?200 
is cosmology-dependent but, in any case, the implicit as- 
sumption Zf = Zobs is a greater source of uncertainty. In 
particular, there is a systematic trend for the discrepancy 
between these two quantities to vary with system size, in 
accordance with a hierarchical structure formation scenario, 
in which the smallest haloes form first. The consequences of 
this effect are addressed in section lTT^ Given the local nature 
of our sample, the assumed cosmology has little effect on our 
results. For example, comparing the values of luminosity dis- 
tance obtained for qo — and qo = 0.5: the difference is less 
than 5% for our most distant cluster {z = 0.208), dropping 
to less than 2% for z < 0.1 (i.e. for 94% of our sample). 

Length scales in the cluster models are defined in a 
cosmology-independent form, with the core radius of the gas 
density expressed in arcminutes and the temperature gradi- 
ent in equation |5| measured in keV per arcminute. The con- 
tributions to the cluster X-ray flux, in the form of discrete 
line emission from highly ionized atomic species in the IGM, 
are handled differently between the different sub-samples. 
However, in all cases the gas metallicity was measured di- 
rectly in the analysis and hence this emission has, in effect, 
been decoupled from the dominant bremsstrahlung compo- 
nent, which we rely on to measure the gas density and tem- 
perature. 

The key advantage of quantifying gas density and tem- 
perature in an analytical form, is the ability to extrapolate 
and interpolate these and derived quantities, like gas frac- 
tion and overdensity, to arbitrary radius. Consequently, the 
virial radius and emission-weighted temperature can be eval- 
uated in an entirely self-consistent fashion, and thus we are 
able to determine the above quantities at fixed fractions of 
R200 , regardless of the data limits. 

Clearly, where this extrapolation is quite large (e.g. 
at J?20o) there is potential for unphysical behaviour in the 
gas temperature, which is not constrained to be isothermal. 
This is particularly true when steep gradients are involved 
(i.e. large values of a in equation |21 or values of 7 very 
different from unity in equation]^. A linear temperature 
parametrization is most susceptible to unphysical behaviour 
as it can extrapolate to negative values within the virial ra- 
dius. To avoid this problem, we have identified those linear 
T(r) models where the temperature within R200 becomes 
negative. In each case the alternative, polytropic tempera- 
ture description was used in preference, where this was not 
already the best-fitting model. 



3.2 Cooling flow correction 

The effects of gas cooling are well kn own to influen ce the X- 
ray emission from clusters of galaxies (lFabianll^94^ . Cooling 
flows may bo present in as many as 70 per cent of clusters 
iPcrcs ct al. 1998) particularly amongst older, relaxed sys- 
tems, where merger-induced mixing of gas is not a significant 
effect. Consequently we expect cooling fiows to be common 
in a sample of this nature, as we discriminate against ob- 
jects with strong X-ray substructure, which is most often 
associated with merger events. It is possible to infer mis- 
leading properties for the intergalactic gas, both spatially 
and spectrally, if the contamination from cooling flows is 
not properly accounted for. Specifically, gas density core 
radii - and, co nsequently, the /3 index in equation (see 
iNeumann fc Arnaud 19991 for example) - can be strongly 
biased, as can the temperature profile, particularly as cen- 
tral cooling regions have the highest X-ray flux. 

In all of the sub-samples the effects of central cooling 
were accounted for in the original analysis using a variety 
of methods, which are described in the appropriate sections 
below. The final cluster models therefore parametrize only 
the 'corrected' gas density and temperature profiles; thus, 
we have extrapolated the gas properties inward over any 
cooling region, as if no cooling were taking place at all. 



3.3 Markevitch sample 

The sub-sample of Markevitch (hereafter 'M sample') was 
compiled from several separate studies and comprises spa- 
tial and spectral X-ray da ta for 27 clusters of galaxies 
(Mar kevitch et aljlT9 98: Ma rkevitclJl 998: Markevitc h et alJ 
I1999I: iMarkevitch fc Vikhlininlll997l : lMarkevitchlll99ll . Of 
these datasets, 22 are included in our final sample, the re- 
maining systems being covered by one of the other sub- 
samples (the factors affecting this choice are described in 
section |1J. 

To measure the spatial distribution of the gas. X-ray 
images of the clusters were fitted with a modified version 
of equation under the assumption of isothermality, equa- 
tion leads to an equivalent expression for the projected 
X-ray surface brightness, S, given by 



S{r) = S{0) 



1 + 



1 -3/3+i 



(6) 



in terms of projected radius, rp as well as the density core 
radius, rc, and index, /3. This is a modified King function or 
isothermal /3-model (Cavalicre & Fusco-Fcrmiano 1976). For 
all but one of the clusters, data from the ROSAT PSPC were 
used for the surface brightness fitting, as this instrument 
provides greatly superior spatial resolution compared to the 
ASCA telescope (for Abell 1650, no PSPC pointed data were 
available and an Einstein IPC image was used instead). 

Although strictly only appropriate for a uniform gas 
temperature distribution, this approach is valid since, for 
the majority of the clusters in this sub-sample, the expo- 
nential cutoff in the emission lies significantly beyond the 
ROSAT bandpass (~0. 2-2.4 keV). Consequently, the X-ray 
emissivity in this energy range is rather insensitive to the gas 
temperature, and therefore scales simply as the square of the 
gas density. These images were also used directly as models 
of the surface brightness distribution in order to determine 
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the relative normalizations between projected emission mea- 
sures in the different regions for which spectra were fitted 
using ASCA data. 

Gas density data for this sub-sample were provided 
in the form of a King profile core radius and /3 index, 
as derived from PSPC data, using equation El However, 
the density normalization was only available in the form 
of a central electron number density f or a small num- 
ber o f clusters: AbeU 1650 & Abell 399 j.Tones fc FormanI 
Il999^ and Abell 3558, Abell 3266 AbeU 2319 & AbeU 119 
ijMohr. Mathieson fc Ev rard"l99fll). In the original Marke- 
vitch analyses, density normaliz ation data for the remainin g 
systems were taken from lVikhlinin. Forman fc Jonej il999tl . 
in the form of values of the radius enclosing a known over- 
density with respect to the average baryon density of the 
Universe at the observed cluster redshift. It was therefore 
necessary, for this work, to convert these values into central 
electron densities, to provide the necessary normalization 
component in the cluster models. 

Radii of overden sity of 2000, R' , were taken from 
IVikhlinin et alJ (^9^ and were combined with the gas den- 
sity core radii, rc, and (3 indices to determine the density 
normalization, p(0), given that 



P(0) 











1 + 


id] 


Jo 







4-Kr dr = 



2000p(zobs), (7) 



where p(zobs) is the mean density of the Universe at the 
observed redshift of the cluster. The integration was per- 
formed iteratively using a generalisation of Simpson's rule 
to a quartic fit, until successive approximations differed by 
less than one part in 10*. 

The fitted gas density and temperature data for the 
M sample were corrected for the effects of central gas cool- 
ing in the original analyses: the cluster models based on 
these data parametrize only the uncontaminated cluster X- 
ray emission. This was achieved by excising a central region 
of the surface brightness data in the original analysis and, 
for the temperature data, by fitting an additional spectral 
component in the central regions (where required), to char- 
acterise the properties of the co oling gas flux. Full deta ils of 
these methods can be f ound in IVikhlinin et alJ (^9^ and 
iMarkevitch et alj Jl998l) . 

Temperature data for all the clusters in this sub-sample 
were provided in the form of a polytropic index and a nor- 
malization evaluated at 2rc (as deflned in equation This 
radius was chosen as it lay within the fltted data region (i.e. 
outside of any excised cooling flow emission) in all cases. 
These fits results are based on the projected temperature 
profile, but have been corrected for the effects of projection. 
To construct cluster models, it was necessary to calculate 
T(0) from these normalization values, by re-arranging equa- 
tion and substituting r = 2rc to give 



r(0) = T(2rc) 



5 + 



1) 



(8) 



These central normalization values were combined with the 
corresponding polytropic indices and density parameters to 
comprise a 3-dimensional description of the gas temperature 
variation. Errors on all parameters were determined directly 
from the confidence regions evaluated in the original analy- 
ses. 



3.4 Finoguenov sample 

The sub-sample of Finoguenov (hereafter 'F sample') 
comprises X-ray data compiled from several sources, incor- 
porating a total of 36 poor clusters and groups of galaxies 
( FinoEucno^^^onmM]|_^99l |Finogueno^^^Ion«^j2000|; 
Fiiioeucii o^^Dwi^^^Pomnanl 12000 : iFinoguenovet^l] 
120011 which were subject to similar analysis. Of the corre- 
sponding fitted results, 24 were used in the final sample, 
with the remainder taken from one of the other sub-samples 
(the factors affecting this choice are described in section |3J|. 
A combination of ROSAT and ASCA SIS instrument data 
was used to determine the spatial and spectral properties 
of the X-ray emission respectively. 

Values for the King profile core radius and index pa- 
rameter were taken from surface density profile fits (using 
equation |HJ to PSPC images of the clusters, with the ex- 
ception of HCG 51 and MKW 9, where no such data were 
available and a ROSAT HRI and Einstein IPC observation 
were used respectively. A central region of the surface bright- 
ness data was excluded for all systems, to avoid the bias to 
Tc and P caused by emission associated with central gas 
cooling. The best-fitting parameters were used to determine 
the 3-dimensional gas density and temperature distribution, 
via an analysis of ASCA SIS annular spectra, by fitting vol- 
ume and luminosity-weighted values in a series of spherical 
shells, allowing for the effects of projection. In this stage of 
the analysis the central cooling region was included and an 
additional spectral component was fitted to the innermost 
bins, allowing this extra emission to be modelled. A regular- 
isation technique was used to stabilise the fit by smoothing 
out large discontinuities between ad jacent bins. Further de- 
tails o f this method can be found in lFinoguenov fc PonmanI 

To generate cluster models for these objects, it was nec- 
essary to infer a central gas density normalization, as well 
as an analytical form for the temperature profile. Density 
normalization was determined by a core index function (see 
equation fit to the data points, using the {3 index and 
core radius values from the PSPC surface brightness fits. 
This was achieved by numerically integrating equation0(as 
described in section 13.311 between the radial bounds of the 
spherical shells used to determine the fit points, weighted 
by to allow for the volume of each integration element. 
The core radii and /3 index were fixed at their previously 
determined values and p(0) was left free to vary. A best 
fit normalization was then found by adjusting p(0) so as 
to minimize the statistic. Confidence regions for p(0) 
were determined from those values which gave an increase 
in of one. Fitting was performed using the MIGRAD 
method in t he minuit minimization library from CERN 
Jjamejll99l) and errors were found with MINOS, from the 
same package. For the core radius and jS index parameters, 
a fixed error of four per cent was assumed, based on an es- 
timate of the uncertainties in the surface brightness fitting 
iFinoeuenov. Reiprich fc BohringerlbOOll) . 

Since the original density points were measured in units 
of proton number density, it was necessary to convert them 
to electron number density for consistency between the clus- 
ter models. It was also necessary to allow for a known effec- 
tive area offset between the ASCA SIS and ROSAT PSPC 
instruments. This adjustment amounts to a factor of 1.2 
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multiplication to convert from proton number densities in- 
ferred using the former, to equivalent values measured with 
the latter. 

An analytical form for the gas temperature profile was 
obtained from a mass-weighted (i.e. density multiplied by 
the integration element volume, using equation fit to the 
3-dimensional data points, excluding the cooling component. 
For 5 of the coolest groups (IC 4296, NGC 3258, NGC 4325, 
NGC 5129 & NGC 6329), the cold component was not suf- 
ficiently separated from the bulk halo contribution and so 
those central bins that were affected were excluded from the 
analytical fit. The best fit temperature values were subse- 
quently found for both a linear and polytropic description, 
again based on the criterion. 

The parametrization which gave the optimum (i.e. low- 
est) fit to the data points was used, except where this gave 
rise to unphysical behaviour in the model; for three systems 
(AbeU 1060, HCG 94 & MKW 4) the linear T(r) model led 
to a negative temperature within -R200, when extrapolated 
beyond the data region; in these cases a polytropic descrip- 
tion was used in preference. 

3.5 Lloyd-Davies & Sanderson samples 

The sub-sample of Lloyd-Davies (hereafter 'L sample') com- 
prises 19 of the 2 clusters and groups of galaxies anal- 
ysed in the study of lLlovd-Davies. Ponman fc Cano ^ (I2OOOI) 
(AbeU 400 was omitted as it is thought to be a line-of-sight 
superposition of two clusters). Of the corresponding fitted 
results, 14 were used in the final sample, with the remain- 
der taken from either the M or F samples (see section |1J. 
ROSAT PSPC data were analysed for all the objects, with 
data from the wider passband ASCA GIS instrument in- 
cluded to permit the analysis of certain hotter clusters. 

To extend the sample to include individual galaxies and 
also to improve the coverage at low temperatures, an addi- 
tional six objects were analysed - four groups and two early- 
type galaxies (this sub-sample is hereafter referred to as the 
'S sa mple'). T he galaxy groups were drawn from the sam- 
ple of iHelsdon fc Ponman (200Q.) and were chosen as being 
fairly relaxed and having high-quality ROSAT PSPC data 
available. Cooler systems, in particular, were favoured, in 
order to increase the number of low mass objects in the 
sample. The extra objects include two early type galaxies; 
an elliptical, NGC 6482 and an SO, NGC 1553. 

Genuinely isolated early-type galaxies are rare objects, 
given the propensity for mass clustering in the Universe. In 
addition, finding a nearby example of such a system, which 
possesses an extended X-ray halo that has been studied in 
sufficient detail to measure Tx(r), severely limits the num- 
ber of potential candidates. Although NGC 1553 lies close 
to an elliptical galaxy of similar size (NGC 1549) there is no 
evidence from the PSPC data of any extended emission not 
associated with either of these objects, which might other- 
wise point to the presence of a significant group X-ray halo 
(see section r3.5. 11 . NGC 6482, by contrast, is a large ellip- 
tical (Lb ~ 6 X lO^^I/B©) which clearly dominates the local 
luminosity function and which is embedded in an extensive 
X-ray halo (~100 kpc). Its properties indicate that this is 
probably a 'fossil' group (see section Ti. 5. 2|l and as such, its 
properties are expected to differ from those of an individual 
galaxy halo. 



The data reduction and analysis for the S sam- 
ple was performed in a similar way to the study of 
iLlovd-Davies et alJ |200Cf) and a detailed description can be 
found there. The method used involves the use of a spectral 
'cube' of data - a series of identical images extracted in con- 
tiguous energy bands - which constitutes a projected view of 
the cluster emission. A three dimensional model of the type 
described in sect ion r3. II can be fitted directly to these data 
in a forward fitting approach jEvles et al.lll99lH . in order 
to 'deproject' the emission. The gas density and tempera- 
ture are evaluated in a series of discrete, spherical shells and 
the X-ray emissi on in each shell is calculated with a m ekal 
hot plasma code (iMewe. Lemen fc van den Oordll98^ . The 
emission is then redshifted and convolved with the detector 
spectral response, before being projected into a cube and 
blurred with the instrument point spread function (PSF). 
The result can be compared directly with the observed data 
and the goodness-o f-fit is quan tified with a maximum like- 
lihood fit statistic ICas^EzS)- The model parameters are 
then iteratively modified, so as to obtain a best fit to the 
data. 

The contributions to the plasma emissivity from highly 
ionized species, in the form of discrete line emission, is han- 
dled by parametrizing the metallicity of the gas with a linear 
ramp (assuming fixed. Solar-like element abundance ratios), 
normalized to the Solar value. However, the poorer spectral 
resolution of the PSPC requires that the metallicity be con- 
strained to be uniform where only ROSAT data were fitted 
(as for all six extra systems in the S sample) . For those clus- 
ters where ASCA GIS data were additionally analysed in 
the L sample (denoted by a '-I-' in the right-most column of 
Table 0, the gradient of the metallicity ramp was left free 
to vary. 

The use of maximum likelihood fitting avoids the need 
to bin up the data to achieve a reasonable approximation to 
Gaussian statistics: a process which would severely degrade 
spatial resolution in the outer regions of the emission, where 
the data are most sparse. The only constraint on spatial bin 
size relates to blurring the cluster model with the PSF; a 
process which is computationally expensive and a strongly 
varying function of the total number of pixels in the data 
cube. Although the Cash statistic provides no absolute mea- 
sure of goodness-of-fit, differences between values obtained 
from the same data set are x^-distributed. This enables con- 
fidence re gions to be evaluated, for determining parameter 
errors (c.f. iLlovd-Davies et alj|200oll . 

For the S sample, two different minimization algorithms 
were employed to optimise the fit to the data. A modified 
Levenberg-Marquardt method ^lBevingtoJl969^ was gener- 
ally used to locate the minimum in the parameter space. 
Although very efficient, this method is only effective in 
the vicinity of a minimum and is not guaranteed to lo- 
cate the global minimum. In several cases this approach was 
unable to optimise the cluster model parameters reliably 
and a simulated annealing min imization algorithm was used 
JCoffe. Ferrier fc Rogerj|l994^ . However, the disadvantage 
of this technique is the computational cost associated with 
the very large number of fit statistic evaluations required: 
once the global minimum was identified, the Levenberg- 
Marquardt method was u sed to determine paramete r errors, 
in an identical fashion to lLlovd-Davies et alj (I2OOOI) . 

In order to determine errors on derived quantities, such 
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as gravitating mass and gas fraction, we adopt the rather 
conservative approach of evaluating the quantity using the 
extreme values permitted within the confidence ranges spec- 
ified by the original fitted parameters. Ifowever, although 
this method tends to slightly overestimate the errors, as can 
be seen from the intrinsic scatter in our derived masses in 
section 17.21 it is not liable to introduce a systematic bias 
into any weighted fitting of these data. 

For those systems in the L and S samples where a cool- 
ing fiow component was fitted, a power law parametrization 
was used to describe the gas temperature and density varia- 
tions within the cooling radius (also a fitted parameter) . To 
avoid unphysical behaviour at i? = 0, these power laws were 
truncated at 10 kpc, well within the spatial resolution of the 
instrument (for NGC 1395 a cut-off of 0.5 kpc was used to 
reflect the much smaller size of its X-ray halo). 

3.5.1 NGC 1553 

The X-ray spectra of elliptical galaxies comprise an emis- 
sion component originating from a population of discrete 
sources within the body of the galaxy, as well as a possi- 
ble component associated with a diffuse halo of gas trapped 
in the potential well. The contributions of these differ- 
ent spectral components vary according to the ratio of 
the X-ray to optical luminosity of the galaxy (Lx/Lb) 
\Kim. Fabbiano fc Trinchierilll992l) . Since we are interested 
only in the X-ray halo of the systems in this work, we favour 
those galaxies with a high Lx/Lb, where the emission can 
be traced beyond the optical extent of the stellar population. 

A 14.5ks PSPC observation was analysed, in which 
the SO galaxy NGC 1553 appears quite far off axis, al- 
though within the 'ring' support structure. Some 2000 
counts were accumulated in the exposure and the emission 
is detectable out to a radius of 4.8 arcmin (21 kpc). Al- 
though its Lx/Lb, of 1.53 x 10""^, does not mark it out as 
a particularly bright galaxy, its X-ray halo is clearly vis- 
ible and uncontaminated by group or cluster emission. In 
fact, this ratio is typical of non-group-dominant galaxies 
(c.f. Hclsdon ot al. 2001j). However, for this reason we ex- 
pect a reaso nable contribution to the X-ray flux from dis- 
crete sources: iBlanton. Sarazin fc IrwinH200j) have recently 
found that diffuse emission only accounted for ~84 per cent 
of the total X-ray luminosity in the range 0.3-1 keV, based 
on a 34ks observation with the ACIS-S detector on board 
the Chandra telescope. 

The PSPC data show evidence of central excess emis- 
sion, which is adequately described by a power law spec- 
trum, blurred by the instrument PSF, with a photon in- 
dex consistent with unity. This was modelled as a separate 
com ponent, so as to decou ple its emission from that of the 
halo. iBlanton et alJ i200 J) flnd evidence of a central, point- 
like source which they flt with an intrinsically absorbed disk 
blackbody model. The spatial properties of the X-ray halo 
are not addressed in their analysis, but in any case the emis- 
sion is only partly visible, due to the small detector area of 
the ACIS-S3 CCD chip. 

3.5.2 NGC 6482 

The elliptical galaxy NGC 6482 is a relatively isolated 
object, which has no companion galaxies more than two 



magnitudes fainter within 1 h'^g Mpc. However, its X-ray 
luminosity is in excess of lO'^^^h^g erg s~^, which is very 
large for a single galaxy. These properties classify this ob- 
ject as a 'fossil' group - the product of the merger of a 
number of smaller galaxies, bound in a common poten- 
tial well JPonman et alJll994t iMulchaev fc Zabludoflll999l : 
IVikhlinin et al.lll999l: Ijones. Ponman fc Forbesll2000l) . Cor- 
respondingly, this system is more closely related to a group- 
sized halo - albeit a very old one (c.f. Ijones et all l2000l) 
- than to that of an individual galaxy. The X-ray over- 
luminous nature of this galaxy {Lx/Lb = 0.048) implies 
that the vast majority of the emission originates from its 
large ( > 100 kpc) halo, with a negligible contribution from 
discrete sources. 

Approximately 1500 counts were accumulated in an 
8.5ks pointing with the PSPC. During the fitting process 
it was found that there was a significant residual feature in 
the centre of the halo, which may indicate the presence of 
an AGN. It was not possible to adequately model this fea- 
ture with either a point-like or extended component and it 
was necessary to excise a central region (radius 1.2 arcmin) 
of the data to obtain a reasonable fit. As a result, the core 
radius was rather poorly constrained and hence was frozen 
at its best-fitting value of 0.2 arcmin for the error calcu- 
lation stage. In addition, the hydrogen column could not 
be constrained and had to be frozen at the galactic value 
(7.89 X 10^° cm~^ ), as determined from the radio data of 
IStarket alJ lll992D . 



4 CONSISTENCY BETWEEN SUB-SAMPLES 

As a consequence of converting the data from the different 
sub-samples into a uniform, analytical format, we are able 
to adopt a coherent approach in our analysis. By extrapolat- 
ing the gas density and temperature profiles, it is possible 
to determine the virial radius and mean temperature (see 
below) self-consistently, and thus independently of the ar- 
bitrary data limits. Of course, this process of extrapolation 
can potentially introduce other biases, and this is discussed 
in section 17.51 below. In some systems, emissivity profiles 
are affected by significant central cooling and we emphasize 
that in our analysis we have eliminated this contaminating 
component in all of our targets, in order to maintain consis- 
tency between the different sub-samples. In this section we 
present the results of an investigation into the consistency of 
our sample and the agreement between the different analysis 
involved. 

Mean temperatures were calculated for each system, 
by averaging their gas temperature profiles within 0.3J?200, 
weighted by emissivity and excluding any cooling flow com- 
ponent (hereafter referred to as Tew; see column 5 in ta- 
ble 0. Fig. shows the temperatures determined in this 
way, from the F & M samples, compared to the correspond- 
ing values taken from the original analyses. The F sample 
(left panel) shows good agreement, although some discrep- 
ancy is expected, due to differences in the prescription for 
obtaining Tow. However, two clusters are clearly anomalous 
- AbeU 2670 and AbeU 2597. The case of A2670 is a known 
discrepancy, arising from an unusually high background in 
the SIS observation. A2597 is an example of the compli- 
cations of a large cooling flow, which is more readily re- 
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solved in the SIS observation than th e GIS data. The values 
of Tew quoted in iFinoeuenov et all i2000l) for these clus- 
ters are a ctually based on PSPC and GIS data respectively 
jHobbs fc Willmore 19971 and [Markcvitch et al. 1991) and 
not on the SIS data analysed in that paper. However, to 
maintain consistency we have used just the SIS data to con- 
struct our model for these clusters. 

The agreement between Tew values for the M sample 
(right panel) is less good, but here differences are to be 
expected: the method used in this work weights the tem- 
perature profile, between 0.3-R200 and zero radius, by the 
emissivity of the gas as determined by extrapolating p(r) 
and T{r) inwards from beyond t he cooling flow region. In 
contrast, iMarkevitch et all il998r) determine a flux weight- 
ing for their mean temperatures based on their estimate of 
the emission measure from the non-cooling gas within the 
core region. For strong cooling flows, this gives a low weight- 
ing to the central values of T{r) compared to those values 
just outside the cooling zone. Since almost all the systems in 
this sample have polytropic indices in excess of one, their gas 
temperatures increase towards the centre, so the differences 
in the spatial weighting give rise to a systematic difference 
between values of Tew determined with the two methods. 
The overall effect of our analysis is actually to correct for the 
consequences of gas cooling, rather than simply to exclude 
the contribution from the cold component to the X-ray flux. 
This amounts to a simple normalization offset - the mean 
of the values of Tew from the M sample is 18 per cent lower 
than that of the values determined in this work. 

To assess the consistency between the different initial 
analyses in our sample, we studied the models derived for 
four clusters which were common to the M, F and L samples 
(AbeU 2199, Abell 496, AbeU 780 & AWM 7), providing a 
direct comparison of methods. Fig.|21shows the temperature 
and density profiles for each of these systems - in each plot 
the different lines correspond to a different analysis result. 
It can be seen that the density profiles show excellent agree- 
ment in all but the very central regions. At the redshift of 
the most distant cluster [z — 0.057, for A780), 1 arcmin cor- 
responds to roughly 60 kpc and hence these differences are 
confined to the innermost parts of the data. Since these are 
all cooling flow clusters, any discrepancies in the core can 
be attributed to differences in the way the cooling emission 
is handled. In any case, the effects of these discrepancies on 
the global cluster properties are small. The temperature pro- 
files show considerably more divergence, and for the clusters 
A780 and AWM 7, the L sample temperature rises with ra- 
dius, in contrast to the M and F sample models . In the case 
of A780, data from a recent Chandra analysis jPavid et alJ 
I2OOII: iMcNamara et al. 200a) indicate that T(r) does in- 
deed show evidence of a rise with radius within the inner 
~200 kpc in the ACIS-S detector data, although the ACIS-I 
temperature profile exhibits a drop in the outer bin, in the 
range 200-300 kpc. 

The discrepancy between the temperature profiles of 
A780 and AWM 7 is exacerbated by the rise with radius 
seen in the L sample models, which has the compounding 
effects of increasing the size of -R200, as well as steepening 
the gravitating mass profile. However, these clusters have 
two of the most extreme rises in T(r) of any system in our 
sample, and only 5 other systems show any significant in- 
crease in temperature with radius. While it is clear that some 



clusters show evidence of a radially increasing temperature 
profile in their central regions, it is unlikely that this will 
continue out to the virial radius. This presents a fundamen- 
tal problem for a monotonic analytical profile, which must 
inevitably find a compromise: in general the fit is driven by 
the central regions, which have a greater flux weighting. In 
the case of A780, the difference in T(r) leads to a factor 
of 3 difference in the total mass within -R200, between the 
models, although this discrepancy is reduced to 60 per cent 
for the mass within 0.37?2oo ■ The corresponding effect on the 
gas fraction is also less severe, since the total gas mass in- 
creases with -R200- However, for A496 - whose temperature 
profile is more typical of the systems in our sample - the 
agreement between the gravitating mass within -R200 for the 
different models is much better, varying by only 40 per cent. 



5 FINAL MODEL SELECTION 

In order to arrive at a single model for each system, we deter- 
mined an order of preference for the sub-samples, to choose 
between analyses, where overlaps occurred. An initial selec- 
tion was made on the basis of unphysical behaviour in the 
models; the linear temperature parametrization is prone to 
extrapolate to negative values within -R200, and so a number 
of models were rejected on these grounds. Of the remaining 
overlaps, we preferentially select those cluster models from 
the L sample, as this represents the direct application of the 
model to the raw X-ray data and hence should be the most 
reliable method. Application of this criterion leaves just four 
remaining systems, where an overlap occurs between the F 
and M samples. These were resolved on an individual ba- 
sis; in each case the analysis of the data which covered the 
largest angular area was chosen. Since the ability to trace 
halo emission out to large radii is critical in this study, this 
amounts to selecting the more reliable analysis. The param- 
eters for each of the final models are listed in tabled 



6 COMPARISON WITH CHANDRA AND 
XMM-NEWTON 

To provide a further cross-check on our results, we present 
here a comparison of our temperature profiles with those 
measured using the recently launched Chandra and XMM- 
Newton satellites. A2199 has been observed with Chandra 
and an analysis o f these data has recently been presented by 
[johnstonc et al] i2002l) . The projected temperature profile 
shows a increasing T(r) from the core out to ~2.2 arcmin 
(78 kpc), where it turns over and fiattens somewhat- albeit 
with only 2 data points. This turnover radius is identical 
to our own "cooling radius" as determined in t he L sa mple 
analysis (see column 11 of Table 0. Johnstone et alJ s de- 
projected T(r) rises continually with radius, but is limited 
to the central ~4 arcmin of the cluster. 

An XM M-Newton observatio n of A496 was recently 
analysed by iTamura et all ll200lll . The projected temper- 
ature profile rises from the core and turns over at roughly 
3.5 (137 kpc) arcmin, in good agre ement with ou r "cool- 
ing radius" of 3.44 arcmin. Although iTamura et alf s depro- 
jected T(r) peaks at a slightly larger radius (of ~5 arcmin), 
it clearly indicates that the temperature drops significantly 
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Figure 1. Comparison of tlie emission- weiglited temperatures from ttiis work, witli those from the original Finoguenov (left panel) and 
Markevitch (right panel) analyses. The dashed line indicates the locus of equality. 



beyond this point, in qualitative agreement with our profile 
in Fig. 121 However, closer comparison with our results in 
the outer regions of the halo is hampered by the fact that 
the data from both the A2199 and A496 observations are 
restricted to the innermost ~8 arcmin. 

In both these cases, the observed emission is dominated 
by flux from the central portion of the halo, where the 
temperature drops towards the core. Previously this phe- 
nomenon was thought to be a cooling flow, although recent 
higher qualit y data have revealed a lack of cool gas in this 
region (see iB5hringer et al.l |2002| . and references therein). 
This component has been either modelled out or excluded 
from our analysis, to allow us to infer the properties of the 
ambient IGM within this region, which accounts for the dis- 
crepancy between the Chandra and XMM-Newton T{r) and 
our profiles in Fig.|5| However, the potential for bias caused 
by any cooling region is limited, since it is confined to a 
small central part of the halo (the median ratio of the cool- 
ing radius to -R200 in our sample is 5 per cent). We note 
that the most distant cluster in our sample (Abell 2163) 
has a sufficiently small angular size (-R200 ~11 arcmin) 
to allow Chandra to be able to observe most of its halo; 
iMarkevitch fc VikhlininI (1200 j) have shown that its temper- 
ature profile agrees reasonably well with the ASCA T{r), 
which was the basis for the model we have used for this 
cluster. 



7 RESULTS 

7.1 Gas distribution 

Fig. 121 shows the variation in the slope of the gas density 
profile with emission-weighted temperature for the sam- 
ple. It can be seen that, for the hottest systems (> 3- 



4 ke V), 13 is consistent wi th the canonical value of 2/3 
(e.g. iJones fc Formanlll984^ . However, below this temper- 
ature the gas profiles become increasing flattened compared 
to self-similar expectation , in agreement with the work of 
iHelsdon fc PonmanI (I2OOOI) . There is a strong correlation be- 
tween f3 and temperature as measured by Kendall's K statis- 
tic, which gives a significance of 5.8a. 

Intriguingly, the galaxy (NGC 1553) and fossil group 
(NGC 6482) - the diamonds in Fig. O - seem to deviate 
from this general trend. Although there are only two points, 
these are the coolest objects in the sample and NGC 1553 
in particular appears to have a value of P more consistent 
with clusters than with groups of a similar temperature. We 
will revisit this issue in the broader context of galaxy scaling 
properties in section I5TT1 

We fitted a straight line, in log space, to the points 
(both including and exclu ding the two galaxies) u sing the 
ODRPACK software package llBoggs et al.ll989lll992l) . to take 
account of parameter errors in both the X and Y direc- 
tions. The dotted line in Fig. O shows the best fit relation 



P = (0.439 ± 0.06)r" 



excluding the galaxies. The 



index is marginally c onsistent with the logarithmic slope 
of 0.26 ± 0.03 found bv lHorner. Mushotzkv fc Scharj (119981) 
for their literature-based sample, spanning the range ~1- 
10 keV. The flatter slope of our data reflects the greater 
number of hotter clusters in our sample, where the relation 
tends to flatten to approximately /3 = 2/3, indicated by the 
dashed line in Fig. 21 The fit also matche s the d ata points 
from the simulations of lMetzler fc EvrardI (^93) , which in- 
clude the effects of galaxy winds on the IGM, albeit with 
their points having a ~25 per cent higher normalization. A 
fit to the entire sample yields a flatter relation, given by 



/3 = (0.482 ± 0.06)T" 



Although the points seem to 



be reasonably well described by a simple power law, there is 
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Figure 2. A comparison of the gas density and temperature profiles in four clusters common to the M (dotted lines), F (solid lines) 
and L (dashed lines) samples. The vertical lines mark the position of -R200 for each of the different models. 
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Figure 3. The gas density slope parameter (/3) as a function of system emission-weigiited temperature. The diamonds represent the 
two galaxies in the sample. The dashed line indicates the canonical value of /3 = 2/3 and the dotted line is the best fit to the points, 
excluding the galaxies. 
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Figure 4. Gas fraction within O.3R200 a function of system temperature. The diamonds represent the two galaxies in the sample. 



a considerable amount of intrinsic scatter in the data - 80 
per cent more than would be expected from the statistical 
errors alone. 

The variation in /3 is also reflected in the gas fraction 
(/gas), evaluated within a characteristic radius of 0.3i?2oo 
(Fig.|lJ. There is a clear trend (significant at the 6cr level, 



excluding the two galaxies) for cooler systems to have a 
smaller mass fraction of X-ray emitting gas. However, the 
galaxy NGC 1553 lies well below the general cluster relation, 
consistent with the coolest groups, apparently at odds with 
its /3 of approximately 2/3. This behaviour is also evident 
in /gas within i?200, shown in Fig. |S] By contrast, the fossil 
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group NGC 6482 exhibits gas fraction properties which are 
consistent with its (3 - i.e. shghtly above groups of a similar 
temperature in both cases. For the whole sample, the be- 
haviour of the gas fraction within -R200 is only slightly dif- 
ferent from that within 0.37i2oo; there still remains a strong 
(5.4a") trend, although there is some evidence of a levelling 
off above ~5 keV, above which the significance of a correla- 
tion drops to 3.2(T. 

This can also be seen in the mean gas fraction within 
-R200 for those systems hotter than 4 keV, which gives 

(0.163±0.01)/i7o^, as compared to /gas = (0.134 iO.Ol)/!^^ 
for the whole sample. Since the errors in the evaluation of 
this quantity are dominated by systematic uncertainties, we 
use an unweighted mean /gas, which is sensitive only to the 
intrinsic scatter in the data. This behaviour suggests that 
virialized objects may not be 'closed systems', in that some 
of their gas might have escaped beyond ii200, particularly 
for the coolest groups. Ifowever, it must be remembered 
that any effects of systematic extrapolation errors could con- 
tribute to the observed trend. 

To understand the behaviour of the gas fraction across 
the sample. Fig. |H| shows how /gas varies with radius, 
grouped into five temperature bins for clarity. Beyond 
~0.2_R2oo, the profiles lie in order of temperature such that, 
at a fixed radius, gas fraction decreases as temperature de- 
creases, mirroring the trend seen in Fig. 2] This is essentially 
a simple normalization offset and demonstrates that the ef- 
fects of energy injection are more pronounced in less massive 
(i.e. cooler) systems, particularly below ~3-4 keV, as seen 
in Fig. |S] The general trend is for gas fraction to rise mono- 
tonically (beyond ~0.03 -R200) with radius from ~0.02 in 
the core to around 18 per cent at 7?200, for the richest clus- 
ters {kT > 8 keV). This behaviour demonstrates that the 
distribution of the IGM is not similar to that of the dark 
matter, even in the largest haloes, bu t is significantly m ore 
extended, as previously reported (e.g. lDavid et al]|l995ll . 



7.2 The M - Tx relation 

Since the emission-weighted temperature reflects the depth 
of the underlying potential well which retains the X-ray 
gas, a tight relation between system mass and temperature 
is expected. It can be shown that, for the case of simple 
self-similar scaling, M oc T^/^ (see iMohr fc EvrardI Il997l 
for example). Observations generally reveal a steeper rela- 
tion, however, consistent with a breaking of self-similarity, 
as found in other scaling relations (e.g. L — Tx). 

ASCA temperature profiles, on which we rely in this 
work, have a relatively large systematic uncertainty be- 
cause of the wide mirror PSF. A comparison of the ASCA 
profiles from one of the subsamples used here (that of 
[Markovitch et al. 1998) with recent Chandra and BeppoSAX 



radii (e.e. iDavid et alJ 


2OOII: iMarkevitch & Vikhlininll200ll 


iNevalainen et al.ll200ll 


De Grandi & Molendill2002h. At the 



same time, ASCA temperatures in the regions immediately 
adjacent to the central bins in the cooling flow clusters ap- 
pear to be system atically too high, although within their 
uncertainties (e.e. |Dayid et al] 1200 ll : lArnaud et al.l l200lt 
iDe Grandi fc Molendil 1200211 . Direct comparison is limited 
to a few clusters at present. 



It is important to correct for the effects of any central 
cooling flow when calculating the characteristic temperature 
of a cluster, but is is not obvious how best to achieve this. 
In our analysis below, we employ three different methods 
extrapolating over, or excluding the central region, and also 
weighting with the gas density rather than emissivity. The 
justification for using these three different prescriptions for 
T is as follows: 

(i) emission-weighted, extrapolating over CF: this at- 
tempts to fully correct for the presence of a CF and provide 
an estimate of T in the absence of cooling. 

(ii) emission-weighted, excising cooling region (radii of ex- 
cision are listed in column 11 of Table ^1: this method of 
calculating T more closely matches the CF-corrected, spec- 
troscopic measurements which have been used frequently in 
previous work. 

(iii) mass-weighted, extrapolating over CF: this method 
gives values of T which are more naturally obtained from 
numerical simulations, and is less sensitive than emission 
weighting to the properties in the dense central core. 

We have applied these methods to derive T within two dif- 
ferent radii: 

(a) O.3-R200: the majority of our systems have X-ray emis- 
sion detectable to at least this radius, which is typical of 
group detection radii. 

(b) i?20o: this represents our nominal virial radius, and 
more closely matches the detection radii of rich clusters. 

We thus have six different methods of calculating T, includ- 
ing our default method of emission-weighting T{r) within 
0.3-R200 (i.e. Tew, described above and listed in Table 0. 

We have combined these temperature data with our 
gravitating mass measurements (within both O.3-R200 and 
-R200, as appropriate) to give a total of six M — Tx relations. 
Strictly speaking, the masses we derive should be scaled by 
a factor of (1 -t- zc)^'^^^, to allow for the change in mean 
density of the Universe with redshift. We have chosen to 
omit this adjustment, since Zf is unknown, and the assump- 
tion Zf — 2obs is prone to systematically bias the results, as 
mentioned previously. We note, however, that incorporating 
this correction actually makes very little differ ence to the 
best-fitting parameters jFinoguenov et al]|200lll . 

The set of M — Tx relations is plotted in Fig. to- 
gether with the best-fitting power law in each case. The 
fitting was performed in log space, using the ODRPACK soft- 
ware package, using symmetrical errors in both axes derived 
from the half- widths of the asymmetric errors on the original 
values. The upper section of table |5| lists the parameters of 
the fit lines, together with the corresponding scatter about 
the relation, normalized to that expected from the statisti- 
cal errors alone. A series of 1000 random realisations of the 
data was generated by scattering each point away from the 
best-fitting line, using the Icr errors in both X and Y direc- 
tions. The intrinsic scatter was measured for the real data 
and for each simulated dataset, by summing in quadrature 
the orthogonal distance of each point from the best-fitting 
line. The real scatter was then normalized to the mean scat- 
ter from all the realisations, to give the numbers quoted 
in column 5 of table |5| In each case, the level of scatter is 
fully consistent with the errors, thus justifying the use of 
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Figure 5. Gas fraction within R200 as a function of system temperature. The diamonds represent the two galaxies in the sample and 
the dotted line shows the unweighted mean of the whole sample. 
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Figure 6. Spatial variation of gas fraction within a given radius (normalized to R20o)i grouped by system temperature. The solid 
line represents the coolest systems (including the two galaxies) (0.3-1.3 keV), increasing in temperature through dashed (1.3—2.9 keV), 
dotted (2.9-4.6 keV), dot-dashed (4.6-8 keV) and finally dot-dot-dot-dashed (8-17 keV). 



a weighted, orthogonal distance regression to determine the 
best fit. 

For the emission-weighted temperature and mass within 
0.37?2oo (method A in table I^J, we find a best-fitting rela- 
tion of log (M/Mq) = (12.80 ± 0.03) + (1.92 ± 0.06) x logT 
for the whole sample. Exclusion of the two galaxies has a 
negligible efi'ect on this result. Excision of the cooling region 



(method B) leaves the normalization of the M — Tx rela- 
tion unchanged and increases the index only marginally, to 
1.94±0.06. The use of mass- weighting to evaluate T (method 
C) yields a best-fit which is consistent with those of meth- 
ods A and B. It is therefore clear that, within 0.3f?200, the 
AI — Tx relation shows a significantly steeper logarithmic 
slope than the self-similar prediction of 3/2. The agreement 



© 2002 RAS, MNRAS OOO.ITH^ 



16 A. J. R. Sanderson et al. 



between the diflerent methods for obtaining T demonstrates 
the robustness of this result. The behaviour of the M — Tx 
relation within R200 is similar but with a somewhat less 
steep slope: all three measurements of T (methods D,E & 
F) are consistent in producing a best-fitting power law in- 
dex of ~1.84. The two emission- weighted methods (D & 
E) have identical normalizations, but the effect of using a 
mass-weighting is to increase this value by ~60 per cent. 
Although yielding a flatter slope compared to the M — Tx 
relation within 0.37i2oo, this is still significantly steeper than 
the self-similar prediction. Since T{r) generally drops with 
radius, and more of the emission arises at large radius in 
cooler systems, we expect T for the latter to drop more as 
we move from O.3- R200 to -R200 , henc e flattening Af — Tx. 

The study bv lSato et alJ (|200(J) of 83 clusters, groups 
and galaxies observed with ASCA found a logarithmic slope 
of 2.04 ± 0.42, using total mass within -R200 together with 
temperatures determined by spectral fitting. This value is 
more consistent with the slope for our data within O.3-R200 
than -R200, which may indicate that averaging T within 
0.3-R200 provides a closer match with spectroscopically 
measured temperatures, since real X-ray data are rarely 
detect able out to i?200- iNevalainen. Markevitch fc Formara 
i200(]|) measure a logarithmic slope of 1.79 ± 0.14, calculat- 
ing total mass within an overdensity of 1000 (Miooo), also 
using spectroscopically-derived temperatures. Their normal- 
ization, of log M = 13.15, is intermediate between our values 
for 0.3-R200 and -R200, as expected f or an overdensity of 1000 . 
The slope of 1.78 ± 0.09 found bv lFinoguenov et alJ (|200ir i 
for a sample of 39 clusters (using M500) is consistent with our 
relation within -R200, and their normalization of 13.28±0.05, 
lies slightly below our own values within this radius. 

While many X-ray studies appear to suggest that the 
slope of the M — Tx relation is steeper than the self-similar 
prediction of 1.5, it has been sugge s ted th at this may an arte- 
fact of the analysis. iHorner et al.l (|^9^ measure a slope of 
1.78 ± 0.05 for a sample of 38 clusters, using the (3 model to 
estimate masses. In contrast, they find a slope of 1.48 ±0.12 
for a smaller sample of 11 clusters, for which they have spa- 
tially resolved temperature profiles. They attribute the dis- 
crepancy to the simplistic assumption of isothermality (see 
section 1731 and confirm the apparently self-similar slope of 
the M — Tx relation with another sample of 27 clusters with 
virial mass estimates. However, the virial mass estimator is 
known to be susceptible to bias from interloper galaxies and 
the presence of substructure. In addition, the X-ray data for 
their 11 cluster sample are taken from the literature, and 
are therefore expected to be correspondingly heterogeneous. 

The difi'erences between the temperatures obtained with 
the different methods can be gauged by studying the right- 
most two columns of table |21 These show the mean and 
standard deviation of the ratios obtained by dividing the 
values of T found with each of methods B to F with the cor- 
responding ones determined using method A. It can be seen 
that the effect of excising the cooling region, as opposed to 
extrapolating over it, results in an average 2 per cent de- 
crease in T, consistent with the general trend for T(r) to 
increase towards the centre. An even larger drop in T- of 8 
per cent - is observed when comparing the mass-weighted 
values (method C) with the baseline set (A), although the 
spread of ratios is increased (cr = 0.14). By averaging over 
the whole of -R200 (D to F), the mean temperature decreases 



compared to method A, by 8 and 9 per cent for methods D 
and E (emission- weighted) , respectively. The similarity be- 
tween these mean ratios reflects the proportionately smaller 
influence of the cooling region-excision when integrating over 
the entire cluster volume. The mass-weighted T within R200 
shows an even greater drop, of 22 per cent (albeit with 
a = 0.41), compared to method A. This is due to gas mass 
dropping off less sharply than luminosity, lending greater 
weight to the outer regions, where the gas temperature is 
generally lower. 

The level of scatter in our M — Tx data is consistent 
with, or smaller than the scatter expected just from sta- 
tistical errors (depending upon the way in which the tem- 
perature is weighted) - i.e. values of ~0. 5-1.0 in column 5 
of Table |5| - suggesting that our error bounds are some- 
what conservative, as previously described (see section HT!^ . 
This conservative approach helps to allow for extra sources 
of error - for example, simulations have shown that de- 
viations from hydrostatic equilibrium introduce a 15-30 
per cent rms uncertainty into hy drostatic mass estimates 
iEvrard. Metzler fc Navarro! ll99Bl . In any case, it can be 
seen that a power law is not an ideal description of the data 
in several of the plots in Fig. |7| This may reflect the domi- 
nance of systematic effects when extrapolating out to large 
radii, or could indicate that the data follow a different func- 
tional form. Careful inspection of Fig. |7| reveals some evi- 
dence for a convex shape in a few cases, suggesting that the 
logarithmic slope steepens gradually from the cluster to the 
group regime. A con vex M — Tx relation was predicted by 
the simulations of iMetzler fc Evrar d (1997), but only where 
the input energy provided by galaxy winds is assumed to be 
fully retained as thermal energy in the IGM: their simula- 
tions do not show this behaviour in practice, as the extra 
energy is predominantly expended in doing work redistribut- 
ing the gas within the potential. 

More recently, IDos Santos fc Dorel i2002l) have devel- 
oped a purely analytical model, which predicts a convex 
M -Tx relation of the form: M = MqT^ [l -I- ^] with 
To = 2 keV. This leads to a curve with a self-similar slope 
of 2/3 at the high mass, smoothly increasing to an asymp- 
totic value of 3 for kT — > 0. Their model includes the effects 
of non-gravitational heating on the pre-virialized IGM, as 
well as shock heating, and is able to reproduce the observed 
L — Tx relation with a similar, curved fit to the data points. 

Comparisons with mass measurements using data from 
the latest X-ray missions are rathe r limited at present. How- 
ever, a recent Chandra study by lAllen. Schmidt fc Fabiaiil 
i2001l) has found M — Tx and L — Tx relations in agree- 
ment with the predictions of self-similarity, albeit from a 
small sample of only six rich clusters. This study is based 
on analysis of both X-ray data and gravitational lensing in- 
formation and finds good agreement between mass estim ates 
derived from the different methods. I Allen et al.l i200l ') find 
a M — Tx logarithmic slope of 1.51 ± 0.27 within a radius 
of overdensity of 2500, which is approximately equivalent to 
0.3-R200- However, their result is not directly comparable to 
our AI — Tx relations A, B & C above, since the overdensity 
profiles in our sample are not self-similar, so -R2500 actually 
corresponds to a different fraction of -R200 for each system. 

To permit a proper comparison, we have also derived 
masses within i?2500 and have fitted these data in an identi- 
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Figure 7. Total mass as a function of temperature for tliree different temperature prescriptions (rows) and measured within 0.3-R200 
(left column) and -R200 (right column). In each case the solid line indicated the best fit power law. See upper half of Table El for further 
details. 
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Table 2. Summary of results for the power law M — Tx fitting using different mean temperature prescriptions and integration radii 
(for measuring both mass and temperature). Primed models in the lower half of the table are isothermal models, which have identical 
mean temperatures to models A-F, but different total masses. The bold row indicates our default M — Tx relation. Notes: "Multiples of 
the statistical scatter expected from the errors alone, 'denotes the mean ratio of kT divided by the corresponding values obtained with 
prescription A and is the standard deviation of these ratios across the sample (these numbers have been omitted from the lower half 
since the ratios depend only on temperature, which is unchanged). 



cal way to our other M — Tx relations. Fig. |H1 shows our re - 
sults, together with the five clusters from I Allen et alj (1200 ll) 
which they use in their M — Tx sample. It can be seen that 
their data points agree we ll with our values at similar tem- 
peratures. Fitting the five I Allen et al.l i200l[) clusters using 
the same regression technique as we employed above, we 
find a best-fitting slope of 1.56 ± 0.16 with an intercept of 
log Af = 13.12 ± 0.15. If the sixth cluster from their sample 
(3C295) is included in the fit, the best-fitting slope increases 
to 1.64 ± 0.15 and int ercept decre ases to 13.04 ± 0.14. This 
cluster was omitted bv lAllen et alJ as it was the only member 
of their sample without a confirmed lensing mass estimate. 
Fitting the clusters from our own samp le which are hotter 
than 5.5 keV for comparison with the lAllen et alJ (1200 ll) 
analysis (their coolest cluster has T = 5.56 keV), we find 
a logarithmic slope of 1.84 ± 0.14, with a normalization of 
log Af = 12 .80±0.13. This is marginally consistent with the 
lAllen et al.l result. 

If the difference in slope between the two samples of hot 
clusters is real, it might be r elated to the dynamical state 
of the samples. The sample of lAllen et al.l includes only the 
most relaxed clusters, where the assumption of hydrostatic 
equilibrium has been independently verified by lensing mass 
estimates. Our own sample is less well controlled, although 
we have excluded objects which are clearly not in equilib- 
rium. On the other han d, it is clear from Fig. |H| that the 
shallower slope from the I Allen et al.l data is a poor match 
to the relation for cooler clusters, whilst the steeper slope of 
1.84 fits rather well across the entire temperature range. 

Previous studies have suggested that the high and low 
mass parts of the whole M — Tx relation may be charac- 
terised by power laws with different slopes. The cross-over 
temperature between th e two regimes is typically ~3 keV 
jFinoeuenov et aLllioOlh . It is not obvious from our data 
that there is such a break in the M — Tx relation, as 
has b een found for the L — Tx relation (e.g . iFairle y et alJ 
l2000l. and references therein) . iFinoeuenov et al.l (l200lfl find 
a steepening of the logarithmic slope, from 1.48 ±0.11 above 
3 keV to 1.87 ± 0.14 below. However, this behaviour may 



simply be a manifestation of a smooth transition with tem- 
perature, masked by a dearth of cool systems in their sam- 
ple, where the steepening slope is most apparent. Mo re high 
quality data of the type presented by lAUen et alj 12001) , 
but covering a wide temperature range, will be required to 
establish whether the Af — Tx is really convex. What our 
results demonstrate clearly, is that either the relation steep- 
ens towards lower mass systems, or its slope is substantially 
steeper than 1.5. 



7.3 The effects of non-isothermality 

To investigate directly the effects of neglecting spatial vari- 
ations in gas temperature, we have generated an additional 
set of isothermal models for our sample, i.e. with a = 0, 
for a linear T(r), or 7 = 1, for a polytropic IGM. We have 
used the values of T already determined for the six different 
methods described above - with associated errors - to de- 
fine the constant value. These isothermal models have then 
been subjected to an identical analysis to the original set, 
in order to provide a fair comparison of results. 

Fig. |5] shows the M — Tx relation for the isothermal 
sample derived using temperatures from method A (referred 
to as A'). It can be seen that the convex shape evident in 
panel A of Fig. |7|is largely absent, and that a tighter rela- 
tion about the best-fitting line is observed. The parameters 
of this power law fit are given in the lower half of Table |21 
together with equivalent data for the other five isothermal 
Af — Tx samples. Within 0.3J?200 the logarithmic slope in- 
creases marginally for the isothermal models, but within the 
errors, for each of the three methods of measuring T. How- 
ever, for the two emission- weighted methods, the normaliza- 
tion increases by ~15 per cent, although it is unchanged for 
the mass-weighted T. Similar behaviour is observed for the 
Af — Tx data evaluated within 7i2oo: the logarithmic slope is 
slightly steepened for the isothermal case, and the normal- 
ization is increased - for the emission-weighted methods - 
by ~30 per cent. However, the mass- weighted normalization 
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Figure 8. Total mass within -R2500 as a function of emission-weighted temperature. The solid line is th e best fitting pow er-law to 
the points above 5.5 keV (dotted vertical line), i.e. excluding the grey points. The diamonds are the data of lAUen et al.l J2001^ and the 
dashed line is our best- fitting power-law to these data. See text for details. 




Figure 9. Total mass as a function of emission- weighted temperature, evaluated within O.3-R20O1 for an isothermal IGM (method 
A'). The solid line represents the best-fitting power law. See lower half of Table 1^ for further details. 
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decreases by 15 per cent, compared to the non-isothermal 
models. 

It is clear from this that the assumption of isothermal- 
ity leads to an overestimate of the total mass within -R200J 
when an emission- weighted meth od is used to calculat e T. A 
similar conclusion was reached bv lHorner et alJ lll999t) . for a 
sample of 12 clusters, who found that isothermality overes- 
timated the mass by a fact or of 1.7 - a result confirmed by 
iNeumann fc ArnaudI il999f) . The latter authors found that 
the cumulative mass within a given radius for an isothermal 
cluster is significantly steeper than that of a cluster with 
a polytropic index of 1.25 (a value typical of the systems 
in our sample - see Table , with the intersection of the 
two occurring at ~0.35i?2oo. Consequently, the isothermal 
assumption over-predicts the mass for 96 per cent of the 
cluster volume. 

Neglecting temperature gradients in the IGM appears 
to have little or no effect on the logarithmic slope of the 
M — Tx relation and, once again, the observed slopes are 
in good agreement between the three different methods 
of calculating T. T his is in contrast to the prediction of 
[llorncr ct al. (Il999l) . who suggested that the assumption of 
isothermality leads to a steepening in the M — Tx slope, 
which would otherwise be self-similar (i.e. 3/2). However, 
they base this conclusion on an analysis of a small sample 
(12 systems), with data drawn from a number of difi'erent 
sources in the literature. We also find that the rms scatter 
about the best fit M — Tx relations is significantly reduced 
in our isothermal models, and fully consistent with that ex- 
pected from the statistical errors. We conclude that a power 
law seems to provide a good description of the M — Tx re- 
lation for an isothermal IGM. 

The overestimation of the total mass for the isothermal 

case leads to a corresponding underestimation in the total 

gas fraction within ii200, shown in Fig. IIUI The unweighted 

_ 3 

mean gas fraction for the whole sample is (0.110±0.01)/i7q^ , 

_ 3 

as compared to (0.134 ± 0.01)/i7Q^ for the non-isothermal 
case. It can also be seen that the scatter about the mean is 
lower for the isothermal case, although the apparent drop at 
~1 keV, seen in Fig. IIUI is still noticeable. The most obvious 
outlier on this graph is the galaxy NGC 1553 (the left-most 
point). For this system, the isothermal model results in a 
significantly lower /gas, which greatly increases its distance 
from the sample mean. 



7.4 Virial radius 

The precise location of the outer boundary of a virialized 
halo is difficult to quantify and is very rarely directly ob- 
servable. The virial radius is dependent on the mean den- 
sity of the Universe whe n the h alo was formed, as well as 
the adopted cosmology iLacev fc Cole 1993). Clearly it is 
important to be able to define this quantity reliably, since 
we assume that self-similar haloes will have identical prop- 
erties when scaled by Rv The radius enclosing a mean over- 
density of 200 (-R200) is proportional to -Rv in any given 
cosm ology - and lies within R y for all reasonable cosmolo- 
gies I Brvan fc NormanI Il99^ - and scales in a identical 
way ( Navarro et alJll995^ . However, previous studies have 
not always been able to determine -R200, and so have re- 
lied on other means to estimate this quantity. A tight rela- 



tionship between Tew and Rv (and hence -R200) is expected, 
as both these quantities refiect the depth of the gravita- 
tional potential well in a virialized halo; self-similarity pre- 
dicts that R^ oa VT ew (c.f. the size-temperature relation, 
iMohr fc Evrardlll997ll . This proportionality has been con- 
firmed in ensembles of simulated clusters, which provide a 
value for the nor malization in the relat ion. One such exam- 
ple is the work of lNavarro et alJ (^9^, who deduce that 

i?2oo = 0.813 ' (1 + z)-h^o^ Mpc. (9) 

However, their simulations only included adiabatic compres- 
sion and shock heating, and did not allow for the effects of 
energy injection. 

The correspondence between our values of R200 as de- 
termined from the overdensity profile (listed in Table0 and 
those calculated with equation |2| is shown in Fig. 1111 It can 
be seen that there is significant deviation from the locus of 
equality between these quantities, marked by the solid line. 
The largest discrepancy is observed in the smallest haloes, 
indicating that the NFW equation significantly over-predicts 
R200 in these systems (the effect of extrapolation bias is ad- 
dressed in section 17.51 . This is to be expected, given that 
Tow for these objects is most likely to be susceptible to bias 
from non-gravitational heating. To explore the reasons for 
the disagreement between the two methods for calculating 
R200, we have examined the role of temperature gradients 
as the source of the scatter, given their importance in cal- 
culating the gravitating mass (see equation |1J. We have de- 
fined a simple, quantitative measure of the departure from 
isothermality, which, as has already been seen, can exert a 
significant influence on scaling properties (section 17.31 . We 
use the ratio T(0)/T(0.3i?20o), as this is very sensitive to the 
presence of a temperature gradient, and the two distances 
involved bracket the region of interest used to calculate Tow . 

The relationship between this quantity and the ratio of 
the measured -R200 divided by the NFW predicted value, is 
shown in Fig. 1121 There is clearly a strong anti-correlation 
between these quantities, significant at the 6.9a level. Even 
with the two most extreme points removed (the left and 
right-most points on the graph), the significance of the rela- 
tion drops only slightly, to 6.4a. It can be seen that the most 
isothermal systems (clustered around the solid vertical line) 
scatter around the line of equality between the two measure- 
ments of -R200. This demonstrates that the NFW formula is 
valid only for nearly isothermal haloes, and that it otherwise 
over-predicts R200 for the most common case of a radially 
decreasing temperature profile. 

T.5 Extrapolation bias 

Our analysis relies on the validity of extrapolating analytical 
profiles fitted to an inner region of the data, in order to com- 
pensate for the emission which is undetected. However, since 
the extrapolation is, in general, greater for smaller systems, 
there is a potential for introducing a systematic bias in our 
fitting. This is particularly true of the slope of the gas den- 
sity profile, which is best constrained by the emission from 
outer regions of the X-ray halo, and which has been found 
to vary significantly with temperature (see section ITTTl . The 
work of IVikhlinin et alJ fl999) has shown that there is ev- 
idence of a slight steepening of the gas density logarithmic 
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Figure 10. Upper panel: Gas fraction within -R200 as a function of emission-weighted temperature (within -R200 — method 'D'). Lower 
panel: Gas fraction within -R200 for an isothermal IGM. The dotted lines show the unweighted mean of the whole sample. 
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Figure 11. Predicted R200 from the NFW formula (equationEJ plotted against measured ii200- The solid line indicates the line of 
equality. 
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Figure 12. The ratio between measured -R200 a^nd R200 from the NFW formula, as a function of the ratio between T(r) at i? = and 
R = 0.3-R200- Error bars have been omitted for clarity. The lines of equality on both axes are marked; the solid line represents the locus 
of isothermality. 



slope with radius. Such behaviour would naturally lead to 
lower values of (3 being inferred for cooler systems, since their 
smaller haloes would be detectable to a smaller fraction of 
J?v, compared to those of more massive clusters. This effect 
could therefore explain part of the observed correlation be- 
tween 13 and Tew To explore the effects of extrapolation on 
our data, we have performed a series of fits to the surface 
brightness profiles of two clusters (Abell 1795 & Abell 2029), 
investigating the role of the outer radius of the fitted data 
in constraining the fit parameters. 

These clusters were selected as they are rich systems 
(and hence relatively unaffected by energy injection), with 
high quality data (~100,000 and 30,000 counts, respec- 
tively), that cover a fairly large angular extent. This al- 
lows us to trace the emission to a large fraction of i?v and 
means we can analyse a spatial subset of the data, with- 
out approaching the resolution limits of the instrument. In 
addition, we have chosen systems which have cooling fiow 
emission confined to as small a region as possible, compared 
to the gas halo core size, so as to minimize the bias this 
contamination can have on our results. 

Since we are aiming in this section only to explore the 
behaviour of the gas density logarithmic slope as a func- 
tion of radius, we have adopted a different approach to that 
described in section 031 Rather than applying a full depro- 
jection analysis to the data, we have fitted 1-dimensional, 
azimuthally-averaged, surface brightness profiles for each 
cluster. This method permits a much more direct investi- 
gation of the gas density index in the outer regions of the 
halo and allows us use a quantitative measure of goodness- 
of-fit, based on the criterion. 

We obtained azimuthally-averaged surface brightness 
profiles for both clusters, from ROSAT PSPC data, in the 
following way. An image of the cluster was extracted in the 



0.2-2.4 keV band, and point sources above 4.5cr significance 
were masked out. Using the master veto rate, the contri- 
bution to the background from particles was subtracted, 
and the image was then 'fiattened' by dividing by the cor- 
responding exposure map, to correct for the effects of vi- 
gnetting. An estimate of the astrophysical background was 
obtained, based on an annulus extracted from beyond the 
cluster emission, with the PSPC support spokes removed. 
Point sources were also masked out from the annulus, and 
the remaining counts were extrapolated across the field and 
subtracted from the source image. Finally, a radial profile 
was extracted - centred on the peak in the X-ray emission 
- in a series of fixed-width annuli, with a minimum of 50 
counts per bin (see Fig. 1131 . 

A King profile function ( equation Ell wa s fitted to the 
data, using the QDP package llTennantl ^l999"l. and all three 
parameters were left free to vary. A small central region of 
the data was excised, to prevent emission from the cooling 
fiow biasing the results. The fitting was repeated for a subset 
of the data, excluding emission beyond a fixed radius, i?outor, 
so as to investigate any systematic variation in /3 with radius. 
The results of these fits are summarised in table |H1 and the 
best-fitting model to the whole image is shown, together 
with the data points, in Fig. 1131 

It can be seen from Table |3 that there is only marginal 
evidence for a systematic trend in /3 with radius for A2029, 
with a large overlap between the confidence regions. How- 
ever, A1795 seems to show a significantly lower j3 (0.70) for 
the innermost 9', as compared to the profile fitted out to 
17'. An explanation of this apparent flattening can be found 
in the core radius, which is significantly smaller - 2.2', as 
compared to 3.7' for the whole profile. This behaviour is an 
artefact of the excised central cooling region, which is com- 
parable in size to the gas core in the IGM. Consequently, 
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is biased when a large part of the outer profile is excluded 
from the fit - and this propagates through to the best-fitting 
/3 value. This is confirmed by the bottom row in Table |21 
which lists the results of re-fitting the innermost 9' with rc 
frozen at the best-fitting value from the full, 17' profile: the 
corresponding /3 (0.83) is, in fact, identical to the slope for 
the whole profile. This problem is particularly pronounced 
for Af795 due to its fairly large cooling flow: the excised 
central region amounts to ~f 5 per cent of -R200- The radius 
of the cooling flow in A2029 is only about fO per cent of 
R200, a value which is more typical of groups, and hence the 
pot ential for bia s in rc is minimal. 

iLewis et al.l i2000l) find a similar tendency for Tc and 
/3 to decrease, when fitting truncated radial profiles of sim- 
ulated clusters. Similarly, their analysis excludes a central 
portion of the data, corresponding to emission from either 
the central cluster galaxy or a cooling fiow. They also claim 
to find a steepening in the logarithmic slope of the gas 
density with radius in the outer regions, as reported by 
IVikhlinin et alJ lll999f) : however, this discrepancy is only ev- 
ident in the vicinity of the virial radius - a region rarely 
probed by observations of even the hottest clusters. In any 
case, it is clear that, even if such a steepening is a significant 
effect, it is not able to introduce large systematic relative bi- 
ases between groups and clusters in our analysis. 



8 DISCUSSION 

It is clear from the scaling properties we have examined, 
that virialized systems do not exhibit self-similar behaviour. 
The (3 — T relation and the gas fraction data reveal a flat- 
tening of the gas density profiles, which is most obvious in 
the group regime. These observations are consistent with 
energy injection into the IGM by non-gravitational means. 
However, three questions remain unanswered. Firstly, what 
caused this heating; secondly, when did it take place and, 
thirdly, is self-similarity broken only below a certain critical 
temperature, or does the transition occur gradually? 

There are three main candidates for the origin 
of the self-similarity br eaking. Both galaxy w inds (e.g. 
i Llovd-Davies et aljl200oll and AGN heating (e.g. IWu et alJ 
I2OOC1) are able to inject energy at roughly the levels required 
to raise the entropy of gas in the central regions of the IGM. 
How ever, the role of gas cool ing in imposing an entropy fioor 
(e.g. iMuanwong et aLll200j) could also be significant. Al- 
though our results in this paper provide no means of dis- 
criminating between the first two options, our data do allow 
us to address the viability of cooling. The cooling hypothesis 
is offered some support by our gas fraction results: the vari- 
ation in /gas with Tow shown in figures |S] & |H1 is consistent 
with the loss of gas required in the cool systems if cooling is 
to have a significant effect. On the other hand, there is no 
evidence that groups have an excessive total Lb compared 
to clusters, which wou ld be expected if the coo ling gas ul- 
timately formed stars iHelsdon fc Ponmanll20o3) . This may 
indicate that the /gas trend is due to gas being displaced to 
larger radius by the effects o f energy injection, as dem on- 
strated by the simulations of iMetzler fc EvrardI (^93), for 
example. 

However, it should be noted that simulations which 
model the effects of cluster mergers can generate a sim- 



ilarly extended gas distribution, with out energy injec- 
tion from non-gravitati onal processes iNavarro fc Whitd 
119931: iPearce et alJll994h . Similarly, while the large intrin- 
sic scatter in the P — T relation may reffect a difference 
in the level of energy injection between haloes, it can 
also be attributed to the effect s of hierarchical assembly: 
ICavaliere. Menci fc TozzH I^Q^) predict a /? — T relation 
with a 2(T scatter envelope, resulting from merging histories, 
which is in qualitative agreement with the data in Fig. |H] 

The second question raises two possibilities; either the 
heating took place prior to halo collapse (so called 'exter- 
nal' or 'pre-'heating), or most energy injection occurred af- 
ter virialization (internal heating). Although a rather sim- 
plistic distinction, these two scenarios will manifest them- 
selves in different ways on cluster properties, given the ef- 
fect of injecting energy into a medium in which a signifi- 
cant density gradient has already been established. For ex- 
ample, preheating will tend to weaken the shock bound- 
ary and move it outwards as c ompared to internal heating 
jTozzi. Scharf fc Normanll2000l) . Since this boundary marks 
the point where the gas fraction fades into the universal 
value, this amounts to an observable signature. The evidence 
from /gas is tentative, given the extrapolation uncertainties, 
but there is indication of a systematic variation in the total 
gas fraction within i?200 with Tew • Observations of the X-ray 
background suggest that the heating phase took place over 
a time-scale of ~10'^ yr l,Pen 1999) . although the epoch of 
energy injection is n ot constrained. Conversely, observations 
of the entropy floor jPonman et alll999l : lLlovd-Davies et alJ 
l2000f) place a strict upper limit on Zpreh, the redshift at which 
the preheating epoch could have taken place, of Zproh ^ 10. 

The question of on what mass scale the effects of self- 
similarity breaking occur is more readily answered with our 
large sample. Whereas previous work has pointed to a sharp 
transition between groups and clusters, our results do not 
offer muc h support for this hypothes is. The analytical ap- 
proach of IDos Santos fc Dord i2002t) shows that the two 
mass regimes can be unifled with a simple model that incor- 
porates energy injection from non-gravitational processes. 
Their predicted scaling relations show a gradual steeping of 
M — Tx and L — Tx, with decreasing temperature, and indi- 
cate that accretion shocks cannot be completely suppressed 
in groups. Our measured M — Tx data offer some support 
for a convex relation, as opposed to a broken power law, but 
the scatter is rather large. Whilst higher quality data from 
XMM-Newton and Chandra will doubtless shed some light 
on this issue, the greatest uncertainty lies in the systematic 
bias associated with extrapolation to -R200 : it is necessary to 
trace the X-ray haloes of nearby groups out to in order to 
resolve the issue satisfactorily, and this calls for observations 
with a wider fleld-of-view. 

8.1 Galaxies vs. groups 

In a CDM Universe, the formation of structure proceeds 
in a bottom-up fashion, with the smallest haloes virial- 
izing initially and subsequent merging activity leading to 
the hierarchical assemb ly of progressively larger haloes (e.g. 
iBlumenthal et al]ll984) . Consequently, the smallest objects 
tend to be older, having collapsed at an earlier epoch. This 
then leads to differences in the scaling properties, as a result 
of the higher density of the Universe at that time (e.g. the 
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Cluster name Temperature fioutcr ^?outor/^?200 


Normalization 


/3 


Tc red. x^/dof 


(keV) (arcmin) 


(ct s"-"- /degree^) 




(arcmin) 



A2029 



A1795 



t 



9.80 



8.54 



rc fixed* 



18 
12 
9 
17 
12 
9 
9 



0.68 
0.45 
0.34 
0.60 
0.42 
0.32 
0.32 



294 
356 
399: 
173 
209 
401 
169 



+9 
-33 
-74 
-53 
130 
77 
+21 
17 
+41 
30 
+364 
-130 
+5 
5 



0.73 
0.70 
0.68 
0.83 
0.79 
0.70 
0.83 



+0.02 
0.01 
+0.02 
0.02 
+0.03 
0.03 
+0.02 
0.02 
+0.03 
0.03 
+0.04 
0.03 
+0.01 
0.01 



2.51 
2.17 
1.99 
3.72 
3.26 



2.18 



+0.20 
0.20 
+0.24 
0.24 
+0.31 
0.32 
+0.25 
0.24 
+0.33 
0.33 
0.53 



0.59 
3.72 



1.25/23 
1.35/13 
1.72/8 
1.04/50 
1.18/30 
1.09/18 
1.49/19 



Table 3. Summary of results for the 1-dimensional surface brightness fitting within different radii. Errors are Icr. A central region of 
the data was excluded to avoid contamination of by cooling flow emission; the radii of exclusion were 2.7' & 4' for A2029 & A1795, 
respectively. fThis fit was noticeably biased by the excised central region of the data and was repeated: * indicates that the core radius 
was frozen at its previously-determined best-fit value to the whole profile (see text for details). 

o 




0.1 



0.2 



Radius (degrees) 



0.05 



0.1 

Radius (degrees) 



0.2 



Figure 13. Azimuthally averaged surface brightness profiles for A1795 (left panel) and A2029 (right panel). The solid line indicates the 
best-fitting model (see text for details). 



M — Tx relation normalization). However, in the context of 
a preheating prescription, invoked to explain the breaking 
of self-similarity in galaxy groups, the timing of this early 
formation epoch is critical; it is possible that these objects 
virialized prior to the preheating phase. This would give rise 
to behaviour more consistent with massive clusters of galax- 
ies, which are sufficiently large as to be insensitive to the 
effects of energy injection. One possible candidate for the 
origin of this preheating is population III stars, whose for- 
mation precedes even that of galaxies, which may also have 
contr ibuted to the enrichment of the IGM (e.g. iLoewensteinl 
I2OOII) . 

As has already been seen, the properties of the galaxy- 
sized haloes in this sample appear to differ from those of 
groups of a similar temperature. Specifically, the gas den- 
sity index (/?) is rather steeper than expected from a simple 
extrapolation of the f3 — T relation for the whole sample. 
In addition, the SO galaxy NGC 1553 lies noticeably below 
the M — Tx relation. This can be understood in terms of a 
large discrepancy between its redshift of formation and ob- 
servation: there is a bias towards observing nearby objects of 
small mass, but these are likely to have formed earliest of all 
sized haloes. For a given halo mass, the virial temperature 
is proportional to (1 + zt) and so, for NGC 1553, a value of 



zt ~ 3-4 would increase Tew, and push it to the right of the 
M — Tx relation, as observed. Alternatively, the discrepancy 
can be attributed to the effects of non-gravitational heating 
on Tew caused, for example, by outflows from within the 
galaxy itself: a value of Tow — 0.3 keV (compared to the 
actual value of 0.5 keV) would bring the point back on the 
best-fitting M — Tx relation for the whole sample. This ex- 
planation is further supported by its Lx/I/b, which is low 
enough for the stellar population in the galaxy to have sig- 
nificantly infiuenced its X-ray halo (e.g. IPelleerini fc Ciottil 
[199a) . 

In fact, the behaviour of NGC 1553 is sufficiently un- 
usual that it points to an alternative formation mechanism 
to that which is usually invoked for hot gas in groups and 
clusters. For example, the halo could have been built from 
supernova- driven winds , originating in its own stellar popu- 
lation (e.g. ICiotti et al] ll991'). rather than from primordial 
material . This explanation is supported by its Lx / Lb , and 
may explain the anomalously steep l3 index, which would 
otherwise be flattened by the effects of energy injection prior 
to collapse. 

In contrast, the position of the elliptical galaxy 
NGC 6482 on the M - Tx relation is fully consistent with 
the best flt line to the whole sample. This is not surprising. 



© 2002 RAS, MNRAS 000.IHO71 



Gas fraction and the M — Tx relation 25 



as it is likely that this object is a fossil group and hence 
will exhibit group-type behaviour. What is certain is that 
this must be a very old system, given the long time-scale for 
orbital decay and merging of the galaxy members. Although 
not as discrepant as NGC 1553, NGC 6482 shows some sug- 
gestion that its /3 index may be high for its temperature. 
This could indicate that it formed sufficiently early to have 
been relatively unaffected by a phase of energy injection, 
that would have occurred prior to the formation of most 
larger objects. There is no doubt that its halo must be pre- 
dominantly primordial - its Lx/Lb is sufficiently high that 
the influence of its stellar population is negligible, in terms 
of contributing to the total gas mass. 

Given an early formation epoch, consistent with hier- 
archical formation, coupled with a correspondingly larger 
mean density, it is particularly important to consider the 
cooling time of the X-ray gas. At a radius of 0.1i?200, 
NGC 1553 has a gas cooling time (tcooi) of ~6 Gyr and 
NGC 6482 has tcooi — 3 Gyr. For comparison, the cool- 
ing times of the next two coolest systems in our sample - 
HCG 68 and NGC 1395 - are 18 Gyr and ~25 Gyr, respec- 
tively. In the case of NGC 1553, this implies that some form 
of heating mechanism must have prevented significant gas 
cooling, if its halo was formed before z = 1 (corresponding 
to a light travel time of roughly 6 Gyr). This is also sup- 
ported by the fact that there appears to be no evidence of 
strong cooling in the core, where the density is even higher. 
Energy injection from galaxy winds could provide the heat- 
ing mechanism necessary to explain this result, as suggested 
above. For NGC 6482, some mechanism is also needed to 
prevent catastrophic gas cooling, which is not observed even 
in the core, given that this fossil group must be a very old 
system. Once again, its high Lx/Lb rules out a significant 
heating contribution from supernova-driven winds. However, 
there is some evidence o f an AGN component in this galaxy 
l)Goudfrooii et aDll994h . which may provide a suitable re- 
heating mechanism to prevent the establishment of a cooling 
flow. 



9 CONCLUSIONS 

We have studied the scaling properties of the X-ray emit- 
ting gas and gravitating mass of a large sample of clusters, 
groups and galaxy-size haloes. The 3-dimensional variations 
in gas density and temperature - corrected for contamina- 
tion from cooling flows - are parametrized analytically, al- 
lowing us determine all derived quantities in a self-consistent 
manner. We have derived virial radii and emission- weighted 
temperatures from these models and are able to extrapolate 
the properties of the data to measure them at flxed fractions 
of Rv We have also analysed an identical set of isothermal 
models, to investigate the effects of neglecting spatial vari- 
ations in the temperature of the IGM. We summarise our 
main findings below. 

(1) Beta varies strongly with temperature, although there 
is evidence that galaxy-sized haloes do not follow this trend. 
We find a best fit power law relation of the form (3 — (0.44 ± 
06)T°'^''*'''''^ 

(2) There is a 6a correlation between /gas within 0.3-R200 
and temperature, consistent with the variation in /3. This 
trend is weakened only slightly (to 5.4cr) by extrapolating 



the gas fraction to -R200 although, above 4 keV, the signifi- 
cance of this correlation drops to 3.2a. The mean /gas within 

_ 3 

-R200 for the systems hotter than 4 keV is (0.163±0.01)/i7o^ > 

_ 3 

compared to (0.134± 0.01)/i7Q^ for the whole sample. Under 
the assumption of isothermality, the scatter between /gas at 

-R200 and Tow is reduced, as is the normalization, giving a 

_ 3 

mean for the whole sample of (0.110 ± 0.01)h^Q^ . 

(3) Observations of the variation in gas fraction as a func- 
tion of radius in our sample reveal a systematic trend in 
gas fraction with temperature in all but the central regions 
(<0.3-R20o)- This is consistent with the observed trend in 

/gas with Tew 

(4) In our study of the M — Tx relation, we employ two 
additional methods of calculating the average system tem- 
perature, one of which excludes the central region, another 
weighting the temperature with gas density rather than 
emissivity. We apply our three different methods within both 
0.3ii2oo and ii200, for both mass and T, to give a total of 
six Ai — Tx relations. We find that the logarithmic slope of 
the relation is steeper within O.3-R200 but that, even within 
R200, it is inconsistent with self-similarity. There is close 
agreement between the measured slopes found for each of the 
three different prescriptions for T. For the emission-weighted 
T, within i?2oo, we find M = 2.34 x 10^^ x T^'^-Si±o.o6) 
Mq. We find that the effect of assuming isothermality on 
the slope is negligible, but the normalization increases by 
15 and 30 per cent for 0. 3^200 and resp ectively (c.f. 
iNeumann fc ArnaudI IT99I iHorner et ai]|l999l) . indicating 
that the total gravitating mass is significantly overestimated 
in our data when temperature gradients are neglected. In 
addition, the scatter in the relation is reduced (and fully 
consistent with the parameter errors) compared to the non- 
isothermal case. The corresponding best fit relation is given 
by M = 3.02 X 10i=* X r(i «9±o Mq. 

(5) The relation between -R200 and T ew, as deduced from 
simulated clusters jNavarro et al.lll995^ deviates systemat- 
ically from the measured values of -R200, as inferred from 
the overdensity profile. We find a strong negative corre- 
lation between the ratio of the NFW predicted -R200 to 
our measured values and a quantitative measure of non- 
isothermality (T(0)/r(0.37?2oo))- We show that only in the 
absence of a temperature gradient do the methods agree. 

(6) We address the issue of systematic bias associated 
with the extrapolation of the X-ray data to ii200, by fit- 
ting azimuthally averaged surface brightness profiles for two 
clusters, within different outer radii. We flnd no evidence for 
a significant variation in jS with cluster radius and conclude 
that the flatter gas density proflles of cooler systems cannot 
be attributed to the generally smaller angular range over 
which data are available for these objects. 

(7) We flnd that the two galaxies in the sample display 
unusual properties. We have selected these objects on the 
basis of a lack of associated group or cluster halo emission, 
which can contaminate the galaxy halo flux. The SO galaxy, 
NGC 1553 has a steep /3 index and falls to the right of 
the main M — Tx relation, indicative of an early formation 
epoch {z{ ~ 3-4) - which causes haloes of a given mass to be 
hotter than those collapsing at later times. It is also possi- 
ble that Tew for this galaxy may have been artiflcially raised, 
probably by supernova-driven outflows from its stellar pop- 
ulation. The elliptical galaxy NGC 6482 also shows a rather 
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steep gas density profile, but otherwise exhibits group-like 
behaviour, consistent with a classification as a 'fossil' group. 
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